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ABSTRACT 



Here we present the steady state numerical solution and two dimensional hydrodynamic cal- 
culations of supergalactic winds, taking into consideration strong radiative cooling. The two 
possible outcomes: quasi-adiabatic and strongly radiative flows are thoroughly discussed, to- 
gether with their implications on the appearance of supergalactic winds both in the X-ray and 
visible line regimes. 



Subject headings: 
galaxies 

Introduction 



Galaxies: superwinds, Starbursts: optical and X-ray emissions, general - starburst 



Massive starbursts, well localized short episodes 
of violent star formation, are among the intrinsic 
characteristics of the evolution of galaxies. They 
are found both at high and intermediate redshifts 
(Dawson et al., 2002; Ajiki et al., 2002; Pettini et 
al., 1998) as well as in galaxies in the local universe 
(Marlowe et al., 1995, Cairos et al., 2001). 

It has usually been assumed that the energy 
deposition occurs within a region of about 100 pc, 
the typical size of a starburst. However recent op- 
tical, radio continuum and IR observations (Ho 
1997; Johnson et al. 2001, Gorjian et al. 2001) 
revealed a number of unusually compact young 
stellar clusters. These overwhelmingly luminous 
concentrations of stars present a typical radius of 
about 3 pc, and a mass that ranges from a few 
times IO'^Mq to IO^Mq. Clearly these units of 
star formation (super-star clusters) are very differ- 
ent from what was usually assumed to be a typical 
starburst, and we shall put a special emphasis on 
the outflows expected from them. 

The large kinetic luminosity produced by vio- 
lent bursts of star formation is well known to dras- 
tically affect the surrounding interstellar medium. 
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generating the giant superbubbles and supershells, 
detected in a large sample of star-forming galaxies 
(Bomans, 2001). In other cases, the energy input 
rate is able to drive their associated shock waves 
to reach the outskirts of their host galaxy, lead- 
ing, in these cases, to the development of a super- 
galactic wind with outflow rates of up to several 
tens of solar masses per year (see Heckman 2001 
for a complete review). One of the main impli- 
cations of such outflows is the clear contamina- 
tion of the intergalactic medium with the metals 
recently processed by the massive starburst. In 
both cases it is usually assumed that at the heart 
of the starburst, within the region that encom- 
passes the recently formed stellar cluster (Rss), 
the matter ejected by strong stellar winds and su- 
pernova explosions is fully thermalized. This gen- 
erates the large central overpressure responsible 
for the mechanical luminosity associated with the 
new starburst. Within the starburst region, it is 
the mean total mechanical energy hsB and mass 
MsB deposition rates that control, together with 
the actual size of the star forming region Rss, 
the properties of the resultant outflow. The total 
mass and energy deposition rates define the cen- 
tral temperature Tsb and thus the sound speed 
csB at the cluster boundary. As shown by Cheva- 
lier & Clegg (1985; hereafter referred to as CC85) 
at the boundary r — Rsb, the thermal and kinetic 
energy flux amount exactly to 9/20 and 1/4 of the 
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There is however a rapid evolution as matter 
streams away from the central starburst. After 
crossing r = Rsb the gas is immediately acceler- 
ated by the steep pressure gradients and rapidly 
reaches its terminal velocity (Voo ~ 2 csb)- This 
is due to a fast conversion of thermal energy, 
into kinetic energy in the resultant wind. In this 
way, as the wind takes distance to the star clus- 
ter boundary, its density, temperature and ther- 
mal pressure will drop as r^^, r"**/^ and r^^"/"*, 
respectively (CC85). Previous analysis by Heck- 
man et al. (1990) showed that thermal pressure 
radial distributions inside starburst galaxies M82 
and NGC 3256 are in a good agreement with the 
CC85 model. However Martin et al. (2002) found 
that X-ray surface brightness and radial tempera- 
ture gradient observed by Chandra in NGC 1569 
are inconsistent with CC85 predictions. 

Note that the wind is exposed to the appear- 
ance of reverse shocks whenever it meets an ob- 
stacle cloud or when its thermal pressure becomes 
lower than that of the surrounding gas, as is the 
case within superbubbles. There, the high pres- 
sure acquired by the swept up ISM becomes larger 
than that of the freely expanding "free wind re- 
gion" (FWR), a situation that rapidly leads to the 
development of a reverse shock, the thermalization 
of the wind kinetic energy and to a much reduced 
size of the FWR. Thus for the FWR to extend 
up to large distances away from the host galaxy, 
the shocks would have had to evolve and displace 
all the ISM, leading to a free path into the in- 
tergalactic medium and to a supergalactic wind 
with properties similar to those derived by CC85. 
Here we study the true physical properties of such 
well developed FWRs, taking into consideration 
strong radiative cooling. Section 2 compares the 
adiabatic solution of CC85 with our steady state 
solution where radiative cooling is considered. In 
section 2 we also develop an easy criterion to de- 
fine in which cases radiative cooling is to become 
dominant. Section 3 displays two dimensional cal- 
culations that use CC85 as the initial condition 



which evolves until a new steady state is reached. 
Section 4 discusses some of the observational con- 
sequences of such outflows. 

2. The steady state solution 

Following CC85 we assume a spherically sym- 
metric wind, unaffected by the gravitational pull 
caused by the central star cluster or its associated 
dark matter component. The equations that gov- 
ern the steady outflow away from the star forming 
region are: 
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where r is the spherical radius, u{r), p{r) and P{r) 
are the wind velocity, density and thermal pres- 
sure, respectively. Q is the cooling rate (Q = n^A, 
where n is the wind number density and A is 
the interstellar cooling function, Raymond et al. 
1976). At the boundary of the star cluster [Rsb) 
we use the solution of CC85 
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where p is the mean mass per wind particle, k 
is the Boltzmann constant, 7 = 5/3 is the ratio 
of specific heats and cgs is the sound speed at r 
= Rsb, and then solve equations (3)-(5) numeri- 
cally. Note that at r = Rgs the wind Mach num- 
ber M^ = 1. Therefore we use the wind velocity 
as an independent variable at the vicinity of the 
sonic point. Our numerical solutions for Q = 0, 
fully reproduce CC85 adiabatic results. 
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2.1. Radiative cooling in supergalactic 
winds 

A first order of magnitude estimate of whether 
or not radiative cooling could affect the thermody- 
namics of supcrwinds, results from a comparison 
of the radiative cooling time scale 

rcoolir) = (11) 
with the characteristic dynamical time scale 

The wind density, which, as it streams away, de- 
creases as (1/r^), in combination with the inter- 
stellar cooling curve, fully define locally the ratio 
of Tcooi/Tdyn- One can infer then that if the ratio of 
the above two quantities becomes smaller than 1, 
radiative cooling sets in, affecting the thermody- 
namical properties of the flow. These simple esti- 
mates indicate that radiative cooling may become 
more efficient the more compact a star cluster is. 
In fact, it is the density of the free wind what re- 
ally matters. For the same total starburst energy, 
the more compact a starburst is, the larger the 
wind density value and thus the larger the cooling 
rate. The temperature profiles for a 10**^ erg s~^ 
starburst with a radius Rsb = 5 pc and differ- 
ent metallicitics of the wind material arc shown in 
Figure 1, for both the adiabatic and the radiative 
cases. Note that for compact starbursts radiative 
cooling sets in within a radius < 10 Rsb- Despite 
the rapid drop in temperature the velocity and 
the density distributions remain almost unaffected 
by radiative cooling. In all cases the thermal en- 
ergy is rapidly transformed into kinetic energy of 
the wind in the close vicinity of the sonic point. 
This leads to a fast gas acceleration that allows the 
wind gas to rapidly approach its terminal velocity 
Voo » 2 csB. At the same time, if radiative cool- 
ing is considered, the temperature profile strongly 
deviates away from the adiabatic solution forcing 
the gas to soon reach temperature values of the 
order of 10* K. 

One can establish an approximate boundary 
that separates the adiabatic from the strongly ra- 
diative cases, from the condition that TcooI 

be- 
comes equal to Tdyn- The corresponding curves for 
different energy input rates, wind terminal veloc- 
ities and wind metallicitics, as a function of Rsb, 



are shown in Figure 2. If the initial wind param- 
eters {LsB and Rsb) intersect below the corre- 
sponding metallicity curve, cooling will be ineffi- 
cient and deviations from the adiabatic solution 
would be negligible. On the other hand, if the 
initial wind parameters intersect above the corre- 
sponding metallicity curve, radiative cooling is ex- 
pected to become important. Note that the larger 
the wind terminal velocity is, the smaller the re- 
gion of the parameter space covered by the cold 
wind solution. 

3. The time dependent solution 

Several two dimensional calculations using 
CC85 as the initial condition adiabatic flows have 
been performed with the Eulerian code described 
by Tenorio-Tagle & Muiioz Tunon (1997, 1998). 
This has been adapted to allow for the continu- 
ous injection in the most central zones, of winds 
with a variety of Lsb and Msb values. The time 
dependent calculations also account for radiative 
cooling and display the evolution from the adi- 
abatic regime to the more realistic final steady 
state solution when radiative cooling is consid- 
ered. There is an excellent agreement between the 
time dependent calculations and the steady state 
solutions displayed in section 2. 

Figure 3 displays isotemperature contours and 
the velocity field within the central 100 pc^ of a 

supergalactic wind driven by an = 10^^ erg 
s"\ csB = 500 km s-i or a Voo = 1000 km s'^, 
emanating from an Rsb = 20 pc. The solution 
for Z = Z0, as expected from Figure 2, is well 
into the strongly radiative regime. Figure 3 shows 
the rapid increase in velocity from its central value 
to the terminal speed ^ 2 csb- The upper panel 
shows the run of temperature for the correspond- 
ing CC85 adiabatic solution, our initial condition, 
while the central panels depict how radiative cool- 
ing proceeds within the supersonic outflow, lead- 
ing to the formation of cool parcels of gas (T ~ 10* 
K) within rapidly expanding shells of cooling gas. 
After a few 10^ years a new steady state solution 
has been reached (lower panel). There the temper- 
ature of the fast moving wind suddenly plummets 
from 10^ K to IC* K at a distance of 37 pc from 
Rsb- Cooling takes place in a catastrophic man- 
ner, a fact that made us search for possible signs of 
fragmentation, however, the sudden pressure gra- 
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dicnts resulting from cooling arc inhibited by the 
strongly diverging flow, leaving both the density 
and the velocity fields almost unperturbed. Sim- 
ilar results were obtained when higher mctallic- 
ity outflows were assumed. The high metallicities 
are indeed expected from the large inflow of new 
metals into the superwind (see Silich ct al., 2001). 
Figure 4 shows the final steady-state solution for 
a wind with Rsb = 5 pc, Lsb = 10^^ erg and 
Voc = 1000 km s~^ when the assumed mctallicity 
of the wind is equal to IZq, 3Zq, 5Zq and IOZq. 
Clearly the impact of this variable is to favor ra- 
diative cooling even closer to the Ksb surface, fact 
that can be also noticed in Figure 1. 

4. The structure of supergalactic winds 

The results from sections 2 and 3 imply that the 
structure of supergalactic winds aught to be re- 
vised, particularly in the case of powerful sources 
with a small value of Rgs in which radiative cool- 
ing reduces the spatial extent of the X-ray emit- 
ting region. In such cases, the only possibility 
for the origin of an extended hot gas component 
arises from shock heated ISM overrun by the cen- 
tral wind. On the other hand, in lower luminos- 
ity and/or widely spread starbursts, the extended 
hot wind region may also contribute significantly 
to the total soft X-ray emission. Here we center 
our attention on very massive and centrally con- 
centrated starbursts, entities that leave no doubt 
that their newly processed elements will sooner or 
later be driven into the IGM. In such cases the 
central wind is affected by radiative cooling and 
the extended structure of such winds will be un- 
detected in the X-ray regime. 

Supergalactic winds present a four zone struc- 
ture: 1) A central starburst region (a source of 
hard and soft X-rays). 2) A soft X-ray emitting 
zone. 3) The line cooling zone. 4) A region of re- 
combined gas, exposed to the UV radiation from 
the central star cluster. 

Figure 5 displays the radial structure of super- 
galactic winds powered by low {Lsb = 10^^ erg 
s~^) and high (Lsb = 10*^ erg s~^) luminosity 
starbursts as a function of Ksb- In all cases wc 
have assumed a free wind region terminal velocity 
equal to 1000 km s~^ and a IZq metal abundance. 
Solid lines mark the distance at which a lO'^'^crg 
superwind acquires a temperature of 5 x lO^K 



and IC^K for an adiabatic solution. Dashed lines 
present the modified location of the two tempera- 
ture boundaries provided by gas cooling. Dotted 
lines mark the position of the two temperature 
boundaries for energetic starbursts (10^'^ erg s~^). 
Note that in both high and low luminosity cases we 
have assumed the same central temperature, and 
thus the same central sound velocity (csb), both 
cases develop the same terminal speed, leading in 
the adiabatic regime to the same temperature ra- 
dial structure. For the low luminosity radiative 
cases (heavy dashed lines) significant departures 
from the adiabatic solution occur for highly con- 
centrated starbursts (say Rsb < 20 pc) bringing 
the zone of rapid radiative cooling (10^ K < T 
< 5 xlO^ K) closer to the starburst surface. On 
the other hand, the radiative solution for the ener- 
getic cases (dotted lines) shows how both temper- 
ature limits here considered move much closer to 
the starbursts, reducing significantly the extent of 
the X-ray and the line cooling zones. The effect is 
noticeable for all Rsb (from Ipc to 100 pc) values 
here considered, a fact that ought to be taken into 
account in full numerical simulations. 

4.1. The impact of photo-ionization and 
the expected Halpha luminosity 

Figure 5 shows that powerful and compact star- 
bursts are to be surrounded by extended, cold and 
fast expanding envelopes. This continuously in- 
jected halo, however, becomes an easy target of 
the UV photons escaping from the central star 
forming regions, and as it presents an R~^ den- 
sity distribution it should be completely ionized 
by the stellar radiation. This is simply because 
the number of recombinations in such a density 
distribution is, as shown by Franco et al. 1990, 
unable to trap the ionization front. As shown in 
Figure 5, the size of the X-ray and line cooling 
zones become strongly affected by radiative cool- 
ing reducing largely their dimensions in favor of a 
closer development of the outer photo-ionized re- 
gion 4. Note that the recombining zone 3 is also 
an additional source of UV photons capable of ion- 
izing the wind outer extended halo. The number 
of recombinations per unit time expected in the 
photo-ionized, outer free wind region is: 
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^=47ri?2„(i?4)T4,(i?4) + 

where the first term in the right hand side refers 
to the number of recombinations in the rapidly 
cooling zone and the second one is related to the 
number of recombinations in the outer cold enve- 
lope (zone 4). Recombinations should lead to Ha 
and Lya emission. One then can estimate the Hq 
luminosity l^Ha = 1-36 x 10~^^iVrec (Leitherer & 
Heckman, 1995) and Lya luminosity l^Lya ~ 8.74 
Lho! (Brocklchurst. 1971). The expected lu- 
minosities calculated for radiative winds powered 
by 10^^ erg s~^ and 10^^ erg s~^ starbursts are 
presented in Figure 6 as a function of Rsb- 

4.2. The X-ray luminosity 

The X-ray luminosity from the hot wind with 
density distribution n(r) = M^/47r/Uau(r)r^ is 

L^ = 47r/ n'K^{T,Z^yAr, (14) 

where Ax{T,Zw) is the hot gas X-ray emissiv- 

ity and Kx is the X-ray emission cutoff radius 
that corresponds to the cutoff temperature Tcut ~ 
5 X lO^K. This is shown in Figure 5. We then 
use a model velocity field u(r) and the cutoff ra- 
dius Kx to calculate a wind soft (O.lkev - 2.4kev) 
X-ray luminosity. The results of the calculations 
for = 10^^ erg s""'^ and = 10^^ erg s~^ 
supcrwinds are shown in Figure 6. We have to 
conclude that the inconsistency between the free 
wind model predictions and observed supcrwind 
X-ray luminosities, mentioned by Strickland & 
Stevens (2000), becomes even larger if one takes 
into account the modifications provided by radia- 
tive cooling. That is, in the standard approach 
the free wind material itself cannot be the main 
source of diffuse superwind X-ray emission. 

5. Conclusions 

We have reanalyzed the physical properties (ve- 
locity, density and temperature) of supcrgalactic 
winds powered by a constant energy input rate, 



taking radiative cooling into account. Three in- 
put parameters control the properties of the wind: 
the total energy input rate Lsb, the characteris- 
tic scale of star formation (Rsb), and a thcrmal- 
ized gas temperature Tsb (or initial wind veloc- 
ity csb). We have revealed two possible outflow 
regimes: a hot wind, which due to its low power 
and low density remains unaffected by cooling, and 
thus as found by CC85 behaves adiabatically. On 
the other hand, winds driven by powerful and/or 
compact starbursts become rapidly dominated by 
catastrophic radiative cooling. 

Supcrwinds driven by compact and powerful 
starbursts (see Figure 2) undergo catastrophic 
cooling close to their sources, and establish a tem- 
perature distribution radically different from those 
obtained in adiabatic calculations. At the same 
time, both velocity and density radial distribu- 
tions remain almost unchanged given the speed 
of the rapidly diverging flow. 

The rapid fall in temperature as a function of 
r reduces the size of the zone radiating in X-rays 
and decreases the siiperwind X-ray luminosity. At 
the same time, cooling brings the 10** K boundary 
closer to the wind center, and promotes the estab- 
lishment of an extended ionized fast moving enve- 
lope. This should show up as a weak and broad 
(~ 1000 km s~^) line emission component at the 
base of the much narrower line caused by the cen- 
tral HII region. 
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Fig. 1. - The intrinsic temperature profile of su- 
perwinds. The dash-dotted hne shows the tem- 
perature distribution derived from the adiabatic 
solution of CC85. Other curves show the effects of 
radiative cooling for different superwind metallic- 
ities. 

Fig. 2. — Adiabatic versus radiative superwinds. 
The lines divide the region defined by the pa- 
rameter space into strongly radiative superwinds 
(above the lines) and adiabatic solutions (below 
the lines), for different values of the terminal ve- 
locity (Voo) and metallicity. 

Fig. 3. Two dimensional superwinds. Isotem- 
perature contours, with a separation AlogT = 0.1 
and the velocity field (longest arrow = 1000 km 
s^^) within the central 100 pc^ of a superwind 
powered by a 10"*^ erg s~^ starburst with a radius 
^SB = 20pc. The evolution starts from (a) the 
adiabatic solution of CC85 and continues at (b) t 
= 1.3 X lO^'yr, (c) 1.6 x lOV and (d) 1.9 x lO^^yr 
when a new steady state solution is reached. In the 
bottom panel the gas temperature falls to 10'*K 
before reaching 40pc from the center. Distance 
between consecutive tick mark = 25pc in all fig- 
ures. 

Fig. 4. — The metallicity effects. The same 
as figure 3 for a 10**^erg s~^ superwind with an 
^SB = 5pc. The panels show the final temper- 
ature and velocity distributions for the radiative 
steady state solution for (a) Z = Zq, (b) SZq, (c) 
SZq and (d) IOZq. 

Fig. 5. The internal structure of free wind re- 
gions. Solid lines mark the distance at which 
a lO'^^erg s^^ superwind, with a velocity Voo = 
1000km s^^, acquires a temperature of 5 x 10'''K 
and lO^'K function of the initial radius Ksb, 
for an adiabatic solution. Dashed lines present the 
modified location of the two temperature bound- 
aries provided by gas cooling (Z = Zq). Dotted 
lines mark the position of the two tcmpcratiire 
boundaries for energetic starbursts (10"'''erg s^^). 

Fig. 6. — Observational appearance of free wind 
regions. The expected X-ray and luminosity 
of superwinds as a function of the starburst power 
and the radius Rsb- 
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